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The transverse energy and charged particle muhiphcity at midrapidity and their ratio are eval- 
uated in a statistical model with the longitudinal and transverse flows for different centrality bins 
1/^ ' at RHIC at ^/snn ~ 130 and 200 GeV. Full description of decays of hadron resonances is applied 

' in these estimations. The predictions of the model at the freeze-out parameters, which were deter- 

, mined independently from measured particle yields and pr spectra, agree qualitatively well with 

CNJ ■ the experimental data. The observed overestimation of the ratio can be explained for more central 

^ ' collisions by the redefinition of dNch/dri\mid- 
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I. INTRODUCTION 
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■ In the previous paper the extensive analysis of two measured global variables, transverse energy {dEr I dr]\mid) 
and charged particle multiplicity {dNch / dr]\m,id) densities at midrapidity, was delivered. The analysis was done in 
I ■ the framework of the single freeze-out statistical model 0,0)01 for the most central collision cases of AGS, SPS and 
RHIC. Now the same method will be applied in the examination of the centrality dependence of the above-mentioned 
variables and their ratio. The main idea of this method is as follows. Thermal and geometric parameters of the 
' model are determined from fits to the particle yield ratios and pr spectra, respectively. Then, with the use of these 
f-^ ^ parameters both densities, dEx/drj and dNch/drj, can be estimated numerically and compared with the data. The 
I first part of this two-step prescription has been already done for the existing midrapidity data for different centrality 
bins at RHIC at ^/SNN — 130 and 200 GeV |^|^. In the present paper the second step will be performed, namely 
' the estimations of dEj'/dr}\mid and dN / drj\mid for these centralities and comparison with the data reported in 
Refs. [Tj, ^, 9, 10] . The main reason for doing it is that the transverse energy and charged particle multiplicity 
measurements are independent of hadron spectroscopy (in particular, no particle identification is necessary), therefore 
they could be used as an additional test of the self-consistency of a statistical model. There is also an additional 
; I ■ pragmatic reason: predictions of variety of theoretical models were confronted with the data in ,lQj , but none of these 
' models was a statistical model. 

The experimentally measured transverse energy is defined as 



L 

Et = J2e, ■ sinO^ , (1) 

1=1 

where 9i is the polar angle, Ei denotes Ei — rriN {mN means the nucleon mass) for baryons, Ei + mN for antibaryons 
and the total energy Ei for all other particles, and the sum is taken over all L emitted particles [lOl| . 

As a statistical model the single freeze-out model is applied (for details see |^). The model succeeded in the accurate 
description of ratios and pr spectra of particles measured at RHIC H, Q . The main postulate of the model is the 
simultaneous occurrence of chemical and thermal freeze-outs, which means that the possible elastic interactions after 
the chemical freeze-out are neglected. The conditions for the freeze-out are expressed by values of two independent 
thermal parameters: T and iib- The strangeness chemical potential ns is determined from the requirement that the 
overall strangeness equals zero. 

The second basic feature of the model is the complete treatment of resonance decays. This means that the final 
distribution of a given particle consists not only of the thermal part but also of contributions from all possible decays 
and cascades. 
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II. FOUNDATIONS OF THE SINGLE-FREEZE-OUT MODEL 

The main assumptions of the model are as follows. A noninteracting gas of stable hadrons and resonances at 
chemical and thermal equilibrium appears at the latter stages of a heavy-ion collision. The gas cools and expands, 
and after reaching the freeze-out point it ceases. The chemical and thermal freeze-outs take place simultaneously. All 
confirmed resonances up to a mass of 2 GeV from the Particle Data Tables ^3i, together with stable hadrons, are 
constituents of the gas. The freeze-out hypersurface is defined by the equation 



— — — — const . (2) 
The four-velocity of an element of the freeze-out hypersurface is proportional to its coordinate 



The following parameterization of the hypersurface is chosen: 

t = r cosh q;|| cosh a^, = r sinhaj^ cos (/), = r sinha^ sint/), = r sinha|| cosha_L, (4) 
where aii is the rapidity of the element, an — tanh~^(r2/t), and controls the transverse radius: 



r: 



2 _ 



rsinha_L. (5) 



The transverse size is restricted by the condition r < Pmax- This means that two new parameters of the model have 
been introduced, i.e. r and Pmax, which are connected with the geometry of the freeze-out hypersurface. 
The invariant distribution of the measured particles of species i has the form 0, Q 



dNi 



(PpT dy 



p''d<J^,fAp-u), (6) 



where da^ is the normal vector on a freeze-out hypersurface, p ■ u = p^u^ , is the four-velocity of a fluid element 
and fi is the final momentum distribution of the particle in question. The final distribution means here that fi is 
the sum of primordial and simple and sequential decay contributions to the particle distribution. The primordial 
part of fi is given by a Bose-Einstein or a Fermi-Dirac distribution at the freeze-out. A decay contribution is a one- 
dimensional or multidimensional integral of the momentum distribution of a decaying resonance (the exact formulae 
are obtained from the elementary kinematics of a many-body decay or the superposition of such decays, for details 
see Q and the Appendix in ij). The resonance is a constituent of the hadron gas and its distribution is also given 
by the Bose-Einstein (Fermi-Dirac) distribution function. Therefore the final distribution fi depends explicitly on T 
and pb- 

With the use of eqs. lO and the invariant distribution takes the following form: 



dN,, 



.3 



d'^PT dy 

where 



T 



J da\\ J sinhaj^(i(sinhQ:x) J dS, p ■ u fi{p ■ u) , (7) 



p ■ u — rriT cosh ay cosh aj_ — pr cos ^ sinh a±_ . (8) 



Note that the above distribution is explicitly boost invariant. 
The pseudorapidity density of particle species i is given by 



dNi /■ 2 dy dN, _ /■ 2 p dN, 



d-q J drj d'^pr dy J d^pr dy 
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Analogously, the transverse energy pseudorapidity density for the same species can be written as 



f ,2 Pt dy dNi f ^ Ei dNi 

= d pt Ei- — — — = d pt Pt -fT 



dr] J ' p drj d^pr dy J Ei d'^pr dy 

For the quantities at midrapidity one has (in the c.m.s., which is the RHIC frame of reference) 



(10) 



drj 



f ,2 PT dNi 

= d PT -j^ — , 

mid J niT d^pT dy 



(11) 



dET,, 

drj 



' f d^pT PT , i = baryon 



= { I d^pT PT , ^ - anUbaryon 



J d'^PT PT 



dNi 
d'^PT dy ' 



i = others 



(12) 



Note that for the older data from RHIC at ^/sJfN — 130 GeV ■8| the antibaryon case is not distinguished and the 
corresponding expression for the transverse energy density is the same as for all other particles. 
The overall charged particle and transverse energy densities can be expressed as 



dN, 



ch 



dr] 



_ dNi 
mid ^ dri 

ieB ' 



(13) 



dET 
drj 



dEj 



, dr] 



(14) 



where A and B {B C A) denote sets of species of finally detected particles, namely the set of charged particles B 
comprises tt"*", tt", , K~ , p and p, whereas A also includes photons, K^, n and n Q. 



III. SUMMARY OF THE FIT PROCEDURE 



Analyses of the particle ratios and pT spectra at various centralities in the framework of the single freeze-out model 
were done for RHIC in [^|^. Here is the brief summary of this approach. It proceeds in two steps. First, thermal 
parameters T and /is are fitted with the use of the experimental ratios of hadron multiplicities at midrapidity. After 
then two next parameters, r and Pmax, are determined from the simultaneous fit to the transverse-momentum spectra 
of TT*, , p and p. The fits are performed with the help of the method. For the fcth measured quantity R^.^^ 
and its theoretical equivalent R'j!^{ai, ai), which depends on I parameters ai, ai, the function is defined as 



,^(,.,....o.)=i: '"'"-^"r-'--'"'^ (15) 

where cr^ is the error of the fcth measurement and n is the total number of data points. The fitted values of parameters 
mean the values at which has a minimum. 

To determine T and hb the fcth measured ratio of hadron multiplicities at midrapidity is put into eq. (|15|l as the 
measured quantity: 



nexp 



dN^P/dy 



dNj^P 



/dy 



(16) 



In the case of a boost-invariant model (as here), the theoretical equivalent R'j!^{T, ]1b) is given by 
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dNl'^/dy 



Idy 



.jiT,HB) 



(17) 



where riiiT^jiB) is the final density of particle species i calculated for a static gas. The last equality follows from 
the assumption that the temperature and chemical potentials are constant on the freeze-out hypersurface. Then the 
volume of the hypersurface factorizes when the invariant distribution ijHJl is integrated formally over all momenta 
to obtain the total multiplicity TV*'* (for details see j^). The final density means here that it collects all decay 
contributions. Thus the final density of particle species i reads 



n,(r, ^iB) = nf~*"'(T, Mb) + ^ ^^(*' (r, ^iB) , (18) 

a 

where ri^J""'"'^'^'^"-^ {T , /is) is the thermal density of the ath particle species at the freeze-out, p(i, a) is the final number 
of particles of species i which can be received from all possible simple or sequential decays of particle a and the 
sum is over all species of resonances included in the hadron gas. The values of T and /is were fixed for both RHIC 
energies (i.e. ^snn = 130 and 200 GeV) with the use of data for the most central collisions 0,0. In the further 
considerations it is assumed that these values are independent of the centrality. This is reasonable since the very 
weak centrality dependence of the particle ratios has been observed so far. Recent analyses done in 0,^3 have just 
confirmed the above-mentioned assumption. 

The second step of this approach is to determine values of r and Pmax ■ Now the fcth measured quantity i?^^^ is the 
value of dN^^^ / {(PpTdy)\„iid for measured particle species i and its transverse momentum pr, whereas the theoretical 
equivalent R]^{t, pmax) is given by the formula iQ. Note that again the function depends on two free parameters, 
now T and Pmax, since the values of T and pB, which were determined early, have been put into eq. 10 . 

The fitted values of parameters of the model are taken from and are gathered in table Since not all bins 

reported in 0, ITsl . IT^ were examined, the lacking values of the geometric parameters have been obtained from the 
linear approximation between the nearest up and down neighbors. This is justified because the geometric parameters, 
when plotted as a function of centrality, show roughly linear dependence jQ]. Also in tableQlthe corresponding number 
of participants {Npart) is given for each bin. If the division into centrality classes is different for the identified charged 
hadron measurements and the dET/drj and dNch/drj measurements, values of Npart are taken from reports on the 
former. Therefore values of Npart from [T3. [l7l | are listed for PHENIX, whereas values of Npart for STAR are taken 
from its transverse energy measurement analysis (in il6j| numbers of participants are not given). 



IV. RESULTS 



The results of numerical estimations of dNch/dr]\mid, eq. divided by the number of participant pairs for 

PHENIX centrality bins tabulated in table |l] are presented in figs.nand[21for ^snn = 130 and 200 GeV respectively. 
Additionally to the straightforward PHENIX measurements of the total charged particle multiplicity, also the data 
from the summing up of the integrated charged hadron yields are depicted in these figures (more precisely, since the 
integrated charged hadron yields are given as rapidity densities, the transformation to pseudo-rapidity should be done, 
which means the division by a factor 1.2 here, see [lal)- This is because fits of the parameters of the model should 
be done to the same pT spectra which are to be integrated to delivered the charged hadron yields. For PHENIX 
at y^SNN = 200 GeV this is not true since fits were done to the preliminary data 15], but the integrated charged 
hadron yields were delivered in the final report p^ . And as one can see from a comparison between figures in [T^ 
and points from the former are slightly above corresponding points from the latter. To estimate the scale of this 
difference, the sum of integrated charged hadron yields at Npart = 114.2 (the point of the biggest discrepancy between 
the model evaluation and the experimental data, see fig.|2Il has been obtained again from digitizing the data depicted 
in preliminary plots in |l5| . This gives 3.31 of charged particles per participant pair, which is exactly 10% above the 
value obtained from the data given in ^.nd the model evaluation is only 8.5% greater than this number. Also 
in ^3 the feeding of p{p) from A(A) decays is excluded, contrary to 15]. Therefore, the model estimates should 
overestimate the corresponding recalculated experimental values for PHENIX at ^/snn = 200 GeV. To diminish this 
effect, integrated p and p yields delivered in [13, were corrected to include back the feeding. The correction was 
done by the division by a factor 0.65, which is rough average of a pt dependent multiplier used by the PHENIX 
Collaboration (see fig. 4 and eq.(5) in [13). Since all just stated reasons for the "systematic" discrepancies affect the 
absolute values of dNch/dri\„rid, they reveal themselves in the normalization factor in the model ( see eq. JT}). For 
PHENIX at y/SNN = 130 GeV the same spectra were used to fit the model parameters and to obtain the integrated 
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TABLE I: Values of thermal and geometric parameters of the model for various centrality bins taken from |^|^. Marked are 
these bins, for which values of geometric parameters have been obtained by the author from linear approximation between the 
nearest neighbours (see the text for the explanation). 



Collision case 


Centrality [%] 




r [fm] 


Pmax [fm] 


PHJ^^NIX at y/sMN = 13U GeV: 


0-5 


348 


8.20 


6.90 


T = 165 MeV, /xs = 41 MeV 


5-15* 


271 


7.49 


6.30 




15-30 


180 


6.30 


5.30 




30-60* 


79 


4.62 


3.91 




60-92 


14 


2.30 


2.0 


PHENIX at Vsjviv = 200 GeV: 


0-5 


351.4 


7.86 


7.15 


T = 165.6 MeV, fiB = 28.5 MeV 


5-10* 


299.0 


7.48 


6.81 




10-15* 


253.9 


7.10 


6.47 




15-20* 


215.3 


6.71 


6.13 




20-30 


166.6 


6.14 


5.62 




30-40 


114.2 


5.73 


4.95 




40-50 


74.4 


4.75 


3.96 




50-60 


45.5 


3.91 


3.12 




60-70 


25.7 


3.67 


2.67 




70-80 


13.4 


3.09 


2.02 


STAR at y/sNN = 200 GeV: 


0-5 


352 


9.74 


7.74 


T = 165.6 MeV, fiB = 28.5 MeV 


5-10 


299 


8.69 


7.18 




10-20 


234 


8.12 


6.44 




20-30 


166 


7.24 


5.57 




30-40 


115 


7.07 


4.63 




40-50 


76 


6.38 


3.91 




50-60 


47 


6.19 


3.25 




60-70* 


27 


5.70 


2.51 




70-80 


14 


5.21 


1.76 



yields [Tj], so the predictions should agree with the recalculated data in principle. As one can notice from figs. ^ and 
13 the above comments concerning both PHENIX measurements are true. The greatest overestimation is in the most 
peripheral bin and in the mid-centrality region of the measurements at ^^snn — 200 GeV. But even at the worst 
point the model prediction is 17% over the corresponding recalculated experimental value. It is not so bad because 
the absolute value of theoretical dNch/dr] is determined by the factor t^. This means that the 5.3% uncertainty in 
r is enough to cause the 17% uncertainty in dNch/dr]. The second reason for the quantitative disagreement is that 
transverse momentum spectra are measured in limited ranges, so very important low-pr regions are not covered by the 
data. To obtain integrated yields some extrapolations below and above the measured ranges are used. In fact these 
extrapolations are only analytical fits, but contributions from regions covered by them account for about 25 — 40% 
of the integrated yields 0]. It might turned out that these extrapolations differ from the thermal distributions 
supplemented by the distributions of products of decays. 

Generally, the agreement of the model predictions with the data is much better for RHIC at ^/snn = 130 GeV. 
For the case of ^/snn = 200 GeV, only the rough qualitative agreement has been reached and the reasons have just 
been explained. It is also worth to stress once more, that the discrepancy between the directly measured dNch/drj 
and dNch/dr] expressed as the sum of the integrated charged hadron yields can be one of these reasons, especially for 
RHIC at y/SNN = 200 GeV (see figEl this effect has already been notified in backup slides of 0]). The discrepancy 
starts at mid-centrality and rises with the centrality. 

The estimates of dEr / dr]\,nid, eq. H14|l . divided by the number of participant pairs are shown in figs. O and 0] for 
y/SNN = 130 and 200 GeV respectively. The quality of the model predictions for dEx/drj measured by PHENIX is 
the same as for dNch/drj. Again, only the qualitative agreement has been reached in the case of ^/snn = 200 GeV 
but the overestimation is higher and equals 30% at most. The main source of the overestimation seems to be the same 
as in the dNch/drj case, namely the uncertainty of fitting the parameter r. The STAR measurements need separate 
comments. The STAR data were taken not at midrapidity {rj — 0) as in the PHENIX case but at 77 = 0.5 on the 
average 0. But the pr spectra used for fits of the geometric parameters were measured at midrapidity, also in the 
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FIG. 1; dNch/drj per pair of participants versus Npart for RHIC at ^/snn = 130 GeV. Dots denote model evaluations, squares 
the newest PHENIX data |10|. triangles the earlier reported PHENIX data and crosses are the recalculated PHENIX data 
from summing up the integrated charged hadron yields delivered in [l^. 
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FIG. 2: dNch/dr] per pair of participants versus Npart for RHIC at ^/sNN ~ 200 GeV. Dots denote model evaluations, squares 
PHENIX data jl^J and crosses are the recalculated PHENIX data from summing up the integrated charged hadron yields 
delivered in |l7|. 
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FIG. 3: dET/drj per pair of participants versus Npart for RHIC at ^/sNN ~ 130 GeV. Dots and open crosses denote model 
evaluations, triangles and squares are PHENIX data [^|10|. Dots and triangles are for the older definition of Et, t.e. the total 
energy Et is taken also for antibaryons in eq. Q. 
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FIG. 4: dET/drj per pair of participants versus Npart for RHIC at y^sjvjv = 200 GeV. Dots and open crosses denote model 
evaluations, squares and stars are data. Dots and squares are for PHENIX (data from ^10^). open crosses and stars for STAR 
(original data from 9] have been rescaled to = 0, see the text for the explanation). 
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FIG. 5: (dEr /d-q) / (dNch/dri) versus Npart for RHIC at ^/sNN ~ 130 GeV and for the older definition of Et, i.e. the total 
energy Ei is taken also for antibaryons in eq. Dots denote model evaluations, squares are the earlier PHENIX data -gj. 
Crosses denote recalculated PHENIX data points, i.e. the sum of integrated charged hadron yields [T^ have been substituted 
for the denominator in the ratio. 



STAR case [l^- Therefore one should expect some "systematic" discrepancy between predictions and the data on the 
whole. To remove this effect the original STAR data Q have been divided by a factor sin (6*1^=0.5) ~ 0.887. As it can 
be seen from fig.01 the predictions and data agree with each other within errors besides the most central point. But 
even there, the discrepancy does not overcome 15%. Of com'se, the second source of the quantitative disagreement 
could be the uncertainty in r (in fact, as it will be seen, it seems to be the main source of the disagreement in all 
discussed cases). 

Fig-El slso shows that both experimental and theoretical values of the transverse energy per participant pair cor- 
responding to the STAR case are ~ 30% greater than the PHENIX ones. As far as the theoretical estimates are 
concerned this is the consequence of the higher values of pt distributions of pions, kaons and antiprotons measured by 
the STAR Collaboration with respect to PHENIX measurements (it can be seen directly from careful comparison of 
spectra given in 0| and in fact such a comparison was done in j^, see fig. 5 there). However, it is difficult to judge 
without doubts what is the reason for the normalization difference between the experimental data. Additionally, it is 
interesting that this difference decreases with the centrality. One of the reasons could be a different fiducial aperture: 
in the PHENIX case measurements were done for | ry |< 0.38 in pseudo-rapidity and = 44.4° in azimuth, whereas 
in the STAR case the pseudo-rapidity range was < 77 < 1 and the azimuthal coverage A<j) — 60°. In both cases the 
correction factor for fiducial acceptance is given by 1/A?7 • 27r/A0, which implicitly assumes the uniform distribution 
of the raw Et data in the covered pseudo-rapidity range and in the azimuthal angle. Of course, it might not be 
exactly true and the final difference in the data could reflect slight deviations from this expected uniformity. Another 
point is a rather rough technique used by the author to rescaled the STAR data to 77 = 0. The division by a factor 
sin (6*1^=0.5) ~ 0.887 contributes almost 13 percentage points to the mentioned ~ 30% increase in the normalization. 
Anyway, the increase in the normalization of the transverse-energy measurements observed in the STAR data with 
respect to the PHENIX data is consistent with the similar increase in the normalization of the measured spectra. 

Values of the ratio {dET / drf) I {dNch/ dr/) as a function of Npart are presented in figs. 15181 As one can see, the position 
of model predictions is very regular and exactly resembles the configuration of the data in each case, the estimates 
are only shifted up about 10% as a whole. This indirectly proves that the earlier discussed disagreement in estimates 
of dNch / drilmid and dEx / drjlmid has its origin in the uncertainty of fitting the parameter r (the normalization factor 
cancels in the ratio). The observed 10% overestimation of the ratio can be explained, at least for more central 
collisions, by the observed discrepancy between the directly measured dNch/drj and dNch/drj expressed as the sum of 
the integrated charged hadron yields. If the original data points are replaced by the recalculated data such that the 
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FIG. 6: {dEr / drf) / {dNch/ drf) versus Npart for RHIC at y^sjvjv = 130 GeV. Dots denote model evaluatfons, triangles are 
PHENIX data jlLO!]. Crosses denote recalculated PHENIX data points, i.e. the sum of integrated charged hadron yields [l^ 
have been substituted for the denominator in the ratio. 
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FIG. 7: (dEr / dr;) / {dNch / drj) versus Npart for RHIC at ^/snn = 200 GeV. Dots denote model evaluations, triangles are 
PHENIX data 10]. Crosses denote recalculated PHENIX data points, i.e. the sum of integrated charged hadron yields 
have been substituted for the denominator in the ratio. 
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FIG. 8: {dEr / drfj / {dNch / drf) versus Npart for RHIC at ^/snn = 200 GeV. Dots denote model evaluatfons, stars are STAR 
data lyl. 



denominators are sums of the integrated charged hadron yields, then much better agreement can be reached for aU 
but peripheral collisions (see figs. l5l7|l . 



V. CONCLUSIONS 



The single freeze-out model has been applied to estimate transverse energy and charged particle multiplicity densities 
for different centrality bins of RHIC measurements at ^Jsnn — 130 and 200 GeV. These two variables are independent 
observables, which means that they are measured independently of identified hadron spectroscopy. Since model fits 
were done to identified hadron data (particle yield ratios and pt spectra) and transverse energy and charged particle 
multiplicity densities are calculable in the single freeze-out model, it was very tempting to check whether their 
estimated values agree with the data. Generally the answer is yes, at least on the qualitative level. As it has just 
turned out, the main source of the quantitative disagreement is the uncertainty in the value of the parameter r. 
The uncertainty strongly influences both densities since their theoretical equivalents contain the normalization factor 
T^. This conclusion is confirmed by the analysis of the transverse energy per charged particle as a function of the 
number of participating pairs. The overestimation of at most 30% obtained for the absolute value of the transverse 
energy density decreases to the overall overestimation of the order of 10% for the ratio {dEr / drj) / {dN ch / drf) . On 
the quantitative level this means that values of only three parameters of the model are confirmed in the present 
analysis, namely T, and the ratio Pmax/T. Since the ratio is directly connected with the maximum transverse-flow 
parameter (for the derivation see ,1]), 



pmax _ Pmax/T {IQ) 

\/l + (fimaxlTY 

this set of parameters is equivalent to T, ps and which are more commonly used in other statistical models 

describing particle production in heavy- ion collisions {e.g. the blast-wave model [li|). 

To summarize, the single freeze-out version of a statistical model fairly well explains the observed centrality depen- 
dence of transverse energy and charged particle multiplicity pseudo-rapidity densities at midrapidity and their ratio. 
Also the dependence on ^/snn of the above-mentioned variables is well recovered . It should be stressed once more, 
that this model very well reproduces the particle ratios and the transverse-momentum spectra measured at RHIC 
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[^|j,|5|,|^. In fact, the description of the identified hadron data was the original motivation for the formulation of the 
model in U 3 . The results presented in this paper confirm, in an independent way, the general conclusion drawn in 
the above-mentioned references, that the single freeze-out model is a very useful (and simple) tool for the estimations 
of fundamental physical parameters like the temperature, chemical potentials or the size of the matter created at the 
final stages of a heavy-ion collision. This supports the idea of the appearance of a thermal system during such a 
collision. 
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